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Abstract — We study the electromagnetic cou- 
pling of the array of Gaussian and multi-moded 
horn antennas to the dual-reflector submillimeter- 
wave telescope on the ESA PLANCK Surveyor de- 
signed for measuring the temperature anisotropies 
and polarization characteristics of the cosmic mi- 
crowave background. In this paper, we present 
the results of our analysis concerning the mea- 
surement of polarization with tilted off-axis dual- 
reflector Gregorian telescope and the propagation 
of multi-moded beams through such a system. 



I. Introduction 

The dual-reflector submillimeter-wave telescope on 
the ESA PLANCK Surveyor is being designed for 
measuring the temperature anisotropies and polar- 
ization characteristics of the cosmic microwave back- 
ground (CMB). 

PLANCK will carry two highly complex focal 
plane instruments, giving the wide frequency cov- 
erage (30GHz to 850GHz) and high sensitivity (of 
AT/T 10 -6 ) necessary to allow for the clean separa- 
tion of the CMB from confusing foreground sources 
of radiation and the determination of the power spec- 
trum of the CMB anisotropies on angular scales down 
to 5 arcminutes. 

Our research is concerned with one of these key 
instruments, the Far-IR High Frequency Instrument 
(HFI), which will cover six frequency bands centered 
at 100, 143,217, 353, 545 and 857GHz (wavelength 
range from 3mm down to 350/xm) [1]. At the tele- 
scope focal plane the HFI consists of an array of 36 
horn antenna structures (Fig. 1) feeding the ultra sen- 
sitive bolometric detectors which will be cryogenically 
cooled to a temperature of lOOmK. 



One objective of our research is the optimization of 
the HFI optical design and the determination of the 
complete set of HFI beam patterns. The challenge of 
the problem is that the telescope is electrically large 
(D/X = 4300 at A = 350/xm) and consists of two 
essentially defocused ellipsoidal reflectors providing a 
very large field of view at the focal plane in order to 
accommodate the two focal plane instruments. 

The latter has a detrimental effect on the quality of 
the imaging introducing coma and astigmatism even 
at the center of the field of view where the highest 
frequency channels are placed. Furthermore the focal 
surface of least confusion has quite an unusual shape 
and it is extremely important to precisely place the 
phase centers of the HFI horn antennas on this sur- 
face. 

Another objective is the characterization of the 
non-conventional polarization properties of the multi- 
beam telescope system. The CMB polarization is ex- 
pected to be at a level of only 10% of the temperature 
anisotropy quadrupole, and the success of the mea- 
surements will depend crucially on the precise knowl- 
edge of the polarization properties of the telescope. 

II. Formulation of the problem 

While the performance of the antenna structures 
can be thoroughly tested in terrestrial conditions, the 
coupling of the HFI with the telescope is, basically, 
optimized through the computer simulations. Such 
simulations are also needed for the data retrieval from 
the raw PLANCK measurements once in orbit. 

The most intricate part of the problem is the com- 
putation of the polarization patterns of the telescope 
beams from the linearly polarized horns and of the 
power patterns of multi-moded beams of extremely 
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Fig.l. Horn positions on the focal plane as seen from the sec- 
ondary mirror (Xre) p and Yjidp ax es in the plane are directed 
upwards and to the left, respectively, while the scan direction 
over the sky is to the right) 

high frequencies (the beams of 12 modes at 54:5GHz 
and of 30 modes at 857GHz). 

Among various simulation techniques, physical op- 
tics (PO) is the most adequate one for the given pur- 
pose. However, conventional implementations of the 
technique [2, 3] do not fit the size of the problem. 
Commercially available packages are also very lim- 
ited in their capacity to rigorously answer this sort 
of questions. For example, even the best commercial 
software requires about one week of computation of 
the main beam of the telescope at the relatively low 
frequency of 143GHz, while all conventional physi- 
cal optics codes collapse at the highest frequencies of 
545 - 857GHz. 

To solve the problem, we developed a special PO 
code [1] that allowed us to overcome the limitations 
of a generic approach for large multi-reflector systems 
and perform typical simulations of the telescope in 
the order of minutes. With recent improvements, it 
requires only 2 minutes for the fairly accurate PO 
simulation of the telescope beam at the frequency of 
lA'iGHz and about 30 minutes for the beam of 30 
modes at 857GHz using a PC Pentium III (500 MHz) 
under the Linux operating system. 

In this paper, we consider the beam of the Gaussian 
horn HFI- 143-1 designed for the separate measuring 
of the fields of two orthogonal linear polarizations, 
'a' and 'b', which correspond to the electric field at 
the beam axis in the sky tilted with respect to the 
local vertical (defined below) by the angle ip a = —45° 
and Vfc = +45°, respectively (the angle ip is mea- 
sured clockwise from the local vertical as seen from 
the telescope). 

We also provide the results of our simulations of 
the defocusing effect of the 30-mode horn HFI-857-1 
from the old layout of the focal plane unit. 

The electric field at the aperture of the Gaussian 
horn HFI-143-1 is linearly polarized and specified by 
the formula 

E(r) = E cxp(-r 2 /wa) e, < r < a (1) 
where r is the radial coordinate, w a = 2.68mm is 
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Fig.2. Power pattern of the telescope beam HFI-143-1 (power 
pattern is the same for both polarizations) 

the beam waist at the aperture, a = 5.0mm is the 
aperture radius and e = e a ,b is the unit polarization 
vector (E(r) = outside the aperture). 

Both the power and phase patterns of such a feed 
coincide perfectly well with the available experimen- 
tal data. The far-field pattern of this particular feed 
satisfies the edge taper requirement of being below 
-25<iB at the angles 6 > 25°. 

H-857-1 is a 30-mode back-to-back conical corru- 
gated horn [4] with the electric field at the aperture 
specified by the set of hybrid modes similar to those 
used in [1] for the HFI-545-1 horn. The total field 
is defined as a non-coherent superposition of all the 
modes of all polarizations. 

Amplitudes of the aperture fields are normalized so 
that each polarized mode has the total power P m = 1. 

III. Polarization of the Gaussian beams 

Fig. 2 shows the power pattern of the telescope 
beam H-143-1 as projected on the plane normal to 
the telescope line-of-sight at the (0, 0) point (a x and 
a y are the horizontal and vertical axes on the plane, 
respectively, measured in degrees). The beam axis 
is at the point a x = —1.365°, a y — 1.197° which is 
defined as the point of maximum power of the beam. 

The beam is well shaped down to —30dB below the 
maximum and can be approximated by a Gaussian 
function (at this level, the pattern does not depend 
on polarization) with a full beam width of W m i n — 8.4 
arcmin and W max = 9.2 arcmin measured at —3dB. 

The polarization of the beam is generally elliptical 
except precisely at the beam axis where it remains 
linear. In order to achieve the required orientation of 
the polarization pattern in the sky, we should orient 
the polarization vector e properly on the horn aper- 
ture [1]. 

For immediate comparison of polarizations mea- 
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Fig. 3. Deviation of the major axis of polarization ellipse from 
local vertical for the telescope beams HFI-143-la and HFI-143- 
lb with (a) ip a = —45° and (b) -i/ij, = +45°, respectively 

sured by different horns when scanning through the 
sky, we should use easily aligned directions in the 
sky as equivalent reference polarization axes for dif- 
ferent beams. Such directions are the meridians in 
the spherical frame of the telescope, with the tele- 
scope spin axis being the pole (the meridians define 
local verticals at various observation points, while the 
parallels are local horizontals that constitute the or- 
thogonal directions). 

Also, we should properly define the reference axis 
for the polarization vectors e of tilted horns. We de- 
fine the reference axis as the direction of e in the horn 
aperture plane that is projected on the vertical axis 
of the focal plane. The orientation of e is specified by 
the angle </> in the horn aperture plane measured from 
this reference in a clockwise direction when looking 
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Fig.4. Magnitude of the minor semi-axis of polarization ellipse 
in the far field of the telescope beam HFI-143-la (both cases 
■ipa = —45° and i/>{, = +45° are rather similar) 

from the horn to the secondary mirror. 

Using these definitions, we have found that the 
beam of the H- 143-1 horn is polarized at its axis at the 
required angles of (a) tp a = —45° and (b) ^ = +45° 
if the horn polarization vector e is specified by the 
angle <p a = —46.91° and fa = +43.09°, respectively. 

In terms of the coordinate frame Ml of the primary 
mirror used as a global frame for the design specifi- 
cation, the unit polarization vectors found above are 
represented by their components as follows 

e a = ( 0.522212, -0.727795, 0.444532 ) 

e b = ( 0.628105, 0.680820, 0.376787 ). (2) 

The deviation of the major axis of the polarization 
ellipse from the local meridian as a function of the 
observation point within the beam in each of the two 
cases is shown in Fig. 3, (a) and (b), respectively. 

When the polarization vector e is properly oriented, 
the cross-polarized component of the far field mea- 
sured along the respective orthogonal direction in the 
sky is minimized. In such a case, the power pattern 
of the cross-polarized component is, basically, deter- 
mined by the power pattern of the semi-minor axis of 
the polarization ellipse at each observation point of 
the beam, Fig. 4. 

The latter is very much the same for any orienta- 
tion of the polarization vector and can be approxi- 
mated as follows 

Per = PcrO [(0/e Q ) Sin(^ - ^ )] P CXp[(0/# O ) 9 ] (3) 

where P cr o, dB, is the maximum power of the minor 
axis component achieved at the points specified by 
the polar angle 6>o and the azimuthal angles ipo ± 90° 
measured from the center of the pattern which is lo- 
cated at the beam axis (the values of the parameters 
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in this approximation depend on the position of the 
horn considered). 

IV. Simulation of multi-moded horns 

Multi-moded horns are designed for receiving max- 
imum microwave power within the required angular 
resolution consistent with the requirements on the 
beam taper at the primary mirror of a level of — 25dB, 
with an aspect angle of 25 degrees. These are rather 
restrictive and contradictory requirements resulting 
in the big aperture area of the corrugated conical 
horns specially optimized for the given application [4] . 

Since the phase front of the aperture field of a coni- 
cal horn is convex, the effective focal center is located 
inside the horn at a distance Rc from the aperture. 
The horn position is specified by the aperture refocus 
parameter Ra — Rf + Rc where Ra is the distance 
from the focal plane to the horn aperture measured 
along the horn axis and Rf is the similar distance 
to the geometrical focus on this axis as found by the 
ray-tracing software. 

One may try to make an estimate of Rc by con- 
sidering the basic mode of the horn as a Gaussian 
beam propagating through the horn aperture. For 
the horn HFI-857-1 (old layout) with the slant length 
L = 28mm and the aperture radius a = 2.55mm, the 
estimate is Rc = 11.9mm while Rf = 9.1mm that 
results in Ra = 21.0mm. 

Accurate simulation shows, however, that such an 
estimate is very misleading. When computing the 
gain G defined as G = 101og(47rP(0)/P ) where P(0) 
is the power of the total multi-moded field on the 
beam axis and Pq is the radiated power of a single 
polarized mode, one can find that G as a function of 
the aperture refocus parameter Ra has a maximum 
at Ra = 13.0mm (Fig. 5). 

The accurate value of Ra is different from the 
rough estimate above by ARa = 8mm = 23A that 
results in an extra gain of about 2.5dB. Thus, the 
multi-moded horns should be placed slightly further 
from the secondary mirror compared to the positions 
found by simple estimates. 

V. Conclusions 

A fast physical optics simulator has been developed 
for the analysis of the dual-reflector submillimeter- 
wave telescope on the ESA PLANCK Surveyor. The 
code overcomes the limitations of a generic approach 
for large multi-reflector quasi-optical systems and can 
perform typical simulations of the telescope in the 
order of minutes. 

A study of the power patterns, polarization char- 
acteristics, defocusing effects and modal structure of 
the telescope beams from both Gaussian and multi- 
moded horns has been performed. Analysis of the 
beams from the linearly polarized Gaussian horns 
have shown that the far-field of the telescope is, gen- 
erally, elliptically polarized except precisely at the 
beam axis where linear polarization is preserved. 
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Fig. 5. Defocusing effect of the 30-mode conical corrugated horn 
HFI-857-1, old layout (1- Gaussian beam approximating the 
basic mode, 2- Basic mode, 3- Total field as a non-coherent 
superposition of 30 modes) 

When rotating the polarization vector of the horn 
field about the horn axis, the major axes of the po- 
larization ellipses at the central part of the telescope 
beam rotate by very much the same angle about the 
direction of observation while the power patterns of 
both the co- and cross-polarized components of the 
far field remain basically unchanged. Systematic de- 
viations from this basic rule which occur mainly in 
the peripheral part of the beam can be successfully 
simulated and taken into account. 

Even for the most tilted Gaussian horns located 
at the edge of the horn array, the power associated 
with the minor axes of the polarization ellipses in 
the telescope beam remains at the level of more than 
30dB below the maximum power of the beam. 

In order to achieve maximum angular resolution 
with multi-moded horns, the latter should be located 
slightly further from the secondary mirror compared 
to the estimates based on the Gaussian approxima- 
tion of the basic mode applied to the optimal focal 
positions as found by the ray tracing techniques. 
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